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Abstract 

The singular manifold method from the Painleve analysis can be used to investigate 
many important integrable properties for the nonlinear partial differential equations. 
In this paper, the two-singular- manifold method is applied to the (2+l)-dimensional 
Gardner equation with two Painleve expansion branches to determine the Hirota bi- 
linear form, Backlund transformation. Lax pairs and Darboux transformation. Based 
on the obtained Lax pairs, the binary Darboux transformation is constructed and the 
N X N Grammian solution is also derived by performing the iterative algorithm 
times with symbolic computation. 
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1. Introduction 

Arising from the Painleve analysis proposed by Weiss, Tabor and Carnevale [l| , the singu- 
lar manifold method (SMM) has been successfully used to investigate the typical integrable 
properties for many integrable nonlinear partial differential equations (NPDEs), such as the 
Lax pair [2, l3|, auto-Backlund transformation , nonclassical Lie symmetry [5[ and 

Hirota bilinear formulation [(>!]. In Refs. 0, 0], it is shown the SMM has turned out to be 
capable of obtaining some special classes of solutions for non-integrable NPDEs. However, 
due to the existence of several Painleve expansion branches for some given NPDEs like the 
modified Korteweg-de Vries (mKdV) equation [sl, Sine-Gordon (SG) equation and mod- 
ified Kadomtsev-Petviashvili (KP) equation j^, in this situation the SMM is not feasible to 
exploit the integrable properties of these equations. Therefore, Refs. [s, 10, [ll| have gen- 
eralized the SMM and developed the two-singular-manifold method to uncover information 
about integrable character. 

Different from the usual expansion, the two- singular-manifold method involves two trun- 
cated Painleve expansions at the constant level term, which contains two different singular 
manifolds at a time. This approach has been applied to the mKdV equation 
equation [sl, classical Boussinesq system 
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11| . Mikhailov-Shabat system 



12|, SG 



|lOl |. generalized 

dispersive long wave equation [q], modified KP equation and so on. With this method, 
not only the auto-Backlund transformation and Lax pair can be obtained, but also the Dar- 
boiix transformation can be constructed in terms of the truncated Painleve expansions in 
both the NPDE and its Lax pair [ol, [isl, [l^. In addition, the relationship relating the singular 
manifolds and Hirota r-function can be precisely established 

Permeation of symbolic computation among various fields of science and engineering re 



markably helps the investigations on the nonlinear partial differential equations (NPDEs) 15, 
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r-i 



17|,ll8|. Symbolic computation has increased the ability of a computer to deal with a 



large amount of complicated and tedious algebraic calculations. 

In this paper, by virtue of the symbolic computation, we will investigate the integrable 
properties for the (2+l)-dimensional Gardner equation 0, [sqI 



u 



yy 



dx + 3 au^ 



0, 



:i.ii 



where a and (3 are two arbitrary constants. When m^, = 0, Eqn. (11. ip reduces to the well- 
known (l+l)-dimensional Gardner equation. For a = 0, Eqn. (11.11) is the KP equation, 
while it is the modified KP equation with /? = 0. Therefore, the (2+l)-dimensional Gardner 
equation could be regarded as a combined KP and modified KP equation. Eqn. (II. ip is 



completely integrable in t 
method Q. Refs. 



20 



le sense that it has been solved by the inverse spectral transform 
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2J] have presented its wide classes of analytical solutions 



including the rational solution, quasi-periodic solution , soliton solution and non-decaying 
real solutions. 



2 



In the following sections, with the help of symbolic computation, we will apply the 
two-singular-manifold method to the (2-|-l)-dimensional Gardner equation to determine the 
Hirota bilinear form, Backlund transformation and Lax pairs. Based on the obtained Lax 
pairs, we will construct the binary Darboux transformation and perform symbolic computa- 
tion on the iterative algorithm to generate the Grammian solutions. 

2. Hirota bilinear form 

To begin with, we rewrite Eqn. (11.11) as the following system 

3 

Ut — Uxxx — Q (3 uUx + -Ci^ Ux — ^Vy + ?) aUxV = (2.1a) 

Vx = Uy. (2.1b) 
Then, we expand the solutions of System (12. ip in a generalized Laurent series 

oo oo 
j=0 j=0 

where x = xi^^U^t), and Uj = Uj{x,y,t), Vj = Vj{x,y,t) are analytical functions in the 
neighborhood of a non-characteristic movable singularity manifold xi^^U^t) = 0? while a 
and b are two integers to be determined. By the analysis of the leading terms, we obtain 

a=l, 6=1, uo = 2e—, Vo = 2e^, (2.3) 

a a 

where e = ±1. It is easy to see that uq and Vq can take two values so that System (12.11) has 
two different Painleve expansion branches. By using two different singular manifolds (p and 
if jo], 3, Q, we take the truncated Painleve expansion at the constant level term 

, 2 f (px ^x\ 1^ . . 

u = u-\ — , (2.4a) 

a \ (f) (f J 

v' = v + ^(^-^), (2.4b) 

where the singular manifold corresponds to e = 1 and to e = — 1, which can also be 
regarded as an auto-Backlund transformation between two different solutions {u' , v') and 
(m, f) for System (12. ip . when singular manifolds and if satisfy the truncation conditions. 
Motivated by Expressions (12. 4p . we introduce the dependent variable transformations 

«=-(--9) =-6og!l , (2.5a) 





^9x 


fx 


a 


V 9 






( 9y 




a 


U 





"V / 

^ ^log| ) , (2.5b) 

y 



to transform System (12.10 into the Hirota bilinear form. Substituting Expressions (12.50 back 
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into System (12. ip . we obtain 

Dt9-f Dlg-f\ ^ ,( Dl9-f \ D.9 ' f , , Dl9-f ( D.9 ' f \ _^fD,g-f^' 



9f 9f Jx \ 9f Jx 9f 9f \ 9f Jx \ 9f 



where D is the well-known Hirota bilinear operator [25 1 

D^D^Dl 9-f = idx- dx^ridy - dy>ndt - d,ygix,y,t)f{x\y',t%,=,,y,=y,,=t. (2.7) 
With symbolic computation, Eqn. (12.61) can be split into 

{Dy + Dl-2(3/aDy) g-f = 0, (2.8) 
(A -Dl + 3 D^Dy -Q/3/aDy) g-f = 0, (2.9) 

which are the Hirota bilinear form of Eqn. (II. ip . By the perturbation technique, one can 
assume the functions / and g in powers of a small parameter e [25| to obtain the multi-soliton 
solutions of Eqn. from Eqns. (ESD and (EJ]). 

It is noted that the key step for the Hirota method is to seek for the suitable dependent 
variable transformation for a given NPDE to be transformed into the Hirota bilinear form. 
If we do not know how to do this, then there is little prospect of being able to use the Hirota 
method. However, the truncated expansion in Painleve analysis can provide us with a useful 
clue in finding such desired transformations. In fact, the Hirota bilinear forms for a large 



class of NPDEs can be obtained in terms of the Painleve truncated expansion 32l, l33|, IS^] 



3. Bilinear Backlund transformation 

For an integrable NPDE, the existence of a Backlund transformation seems to be widely 
accepted 28|. In this section, from Eqns. (12. 8p and (12. 9p . we will derive a bilin- 
ear Backlund transformation between two different solutions u = ^(logg/f)^ and u' = 
^ {logg' / f')^ for Eqn. (11.11) . by considering the following two equations, 

Pl=[{Dy + Dl-2P/aDy)g-f]g'f-gf[{Dy + Dl-2p/aDy)g'-f], (3.1) 
P2= [{Dt-Dl + 3D,Dy-6(3/aDy) g ■ f] g' f 

-g / [ ( A - Dl + 3 D^Dy - 6 (3 /a Dy) g' ■ f] . (3.2) 

With the aid of the Hirota bilinear operator identities (see Appendix flA.ip — flA.Sp ). symbolic 
computation on Eqns. (13. ip and (13.20 yields 

PI = [{Dy + Dl -2P/aDy) g-g']ff'-gg' [{Dy + Dl -2[3/aDy) f ■ f] 

-2D,ign-iD,f-g'), (3.3) 

P2 = 3 D^{D^ g-n- (A f-9')-^ D,{g /') ■ {Dy f ■ g') -3Dy{g /') ■ ( A / ■ 9') 
+ [{Dt-Dl + 3 D,Dy - 6 (3 /a Dy) g ■ g'] f f 

-[{D,-Dl + 3 D^Dy - 6 P/a Dy) f ■ f] g g'. (3.4) 
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Thus, Eqns. (13 .Sp and (13. 4p can be further decoupled into the following equations 



DJ-g' = r]{t)gf', 

v{t)D,g-f' + Dyf-g' + ^{t)gf' = 0, 
[Dy + Dl-2(3/aDy + at)]9-9' = 0, 
[Dy + Dl-2(3/aDy + at)]f-f' = 0, 
[Dt - Dl + 3D,Dy -Q(3/aDy + ({t)] g ■ g' = 0, 
[A -Dl + 3 D,Dy - 6 f3/a Dy + ({t)] / ■ /' = 0, 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 
(3.5e) 
(3.5f) 



where ri{t), 7(t), ^{t) and ({t) are all arbitrary differentiable functions of t. Eqns. (13.51) con- 
stitute the bilinear Backlund transformation for Eqn. (II. ip . from which more complicated 
solutions can be progressively constructed beginning with a seed solution. Additionally, it can 



also be of use for the investigation on other integrable properties [25|, \2M, 1271, l28|, l29|, 
like the nonlinear superposition formula. Lax pair, conservation laws, etc. 
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3l|, 



4. Lax pairs with symbolic computation 

In this section, by the two-singular-manifold method, the Lax pairs of the (2+l)-dimensional 
Gardner equation will be derived. With symbolic computation, we insert Expressions (12. 4p 
into System (12. ip . and get 



(P V (P V 

12 (5 u — 3 + 2vl — Q av — 2wi + QtI + Q aux + ^ fix 
12/9m — 3a^-u^ + 2t'2— 6af — 2w2 + 6r| — Gaw^^ + S V2x 



0, 
0, 



3riy 
3r2y 



V2 V2x 



Wix + 3 Ti Tlx + Vl 

XX 

W2x + ST2 T2x + f 2 XX 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 



with 



A = — (avi — 2j3 + a'^u + aTi), 
2a ^ 

B = - — (^a V2 + 2 P — u — a T2) , 



where Vi, Wi and 



2a 

1,2) are defined as 



(4.6) 
(4.7) 





(pxx 




'^xx 


Vl = 




V2 = 




4>x ' 




fx 








Vt_ 


Wi = 




W2 = 










^x 




(Py_ 


T"2 = 




n = 












^x 



(4.8) 
(4.9) 
(4.10) 
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We use the derivatives of Expression (14.11) with respect to x, y and t, so as to obtain 

A, = A{v2-A-B), (4.11) 

Ay = {AT2), + AB{T2-n), (4.12) 

At = {Aw2), + AB{w2-w{), (4.13) 

B^ = B{vi-A-B), (4.14) 

By = {BT,)^-AB{r,-r,), (4.15) 

Bt = {Bwi)^- AB{W2-Wi). (4.16) 

By virtue of Eqns. (14.61) — ( l4.16p . it is easy to check that the following two relationships are 
satisfied 

{AB)^ = AB{T2-n), (4.17) 

[AB{2 A-2B + vr-V2-^ri-^ rs)], = AB{wi - W2). (4.18) 

Therefore, Eqns. (14.111) — (14. 16p can be simplified as 

A^ = A{v2~ A~ B), (4.19) 

Ay = [AiT2 + B)l , (4.20) 

At = [Aw2 - AB{2 A - 2 B + V, - V2 - 3r, - 3r2)l , (4.21) 

B, = B{vi-A-B), (4.22) 

By = [B{n-A)l, (4.23) 

Bt = [Bwi + AB{2 A-2B + vi-V2-3n-3 ts)]^ . (4.24) 

Through introducing the changes as 

Eqns. (14. 19p — (14.241) can be linearized into 

fc = C(^i-^), (4.26) 

% =C(n-A), (4.27) 

=i^-[wi + A{2A-2B + vi-V2-?>n-?,T2)i (4.28) 

i^t. = i^t{v2-B), (4.29) 

=V^^(r2 + 5), (4.30) 

^+ =V^+[w2-5(2A-2B + t;i-t;2-3ri-3r2)]. (4.31) 

Symbolic computation on Eqns. (14.260 — (14. 3ip with the substitution of Eqns. (14.20 and (14.30 
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gives rise to 

'^y = -^xx - ( «M - — ) V^, (4.32a) 



a 

= 4 + ( 6 a M I il)'^^ + I — ^ Q(3u+-a'^u^ — 'iav + 'iaUx ] il^^ ■, (4.32b) 



and 



V'^ = - ( " ^ - — ) V'^ , (4.33a) 



12/3^^,+ , ^12/32 3 2 2 



W = 4V^j"^^ - ( 6aM ^ j + ( 6/3m + - a u - 3 at; - 3 aw^ ) (4.33b) 



By direct calculation, it is found that Eqn. f 1 1.1 1) can be derived from the compatibility 
conditions = iljfy. Thus, Eqns. (14.321) and (14.331) are two different types of Lax pairs 
of the (2+l)-dimensional Gardner equation. It is noted that through the following gauge 
transformations 

ya a'^ ] 
= exp <^ X -y -t \ r^. 

Lax pairs (14.321) and (I4.33P can be respectively transformed into 

= -r~^ -awT- -/^mF", (4.34a) 
= 4 + 6 a M F~^ + ^6/3M+-a;^M^ — 3af + 3Q; F~ 

+ (3/3m^ + -a/3u2 -3/3t; ) F", (4.34b) 



and 



2 



r+ = F+ -aMF+ + /3uF+, (4.35a) 



3 



F^ = 4 FJt^^ — 6 a u F^^ + |^6/3M+-a^M^ — 3at; — 3a J F^ 

+ (^3/3m^. - ^a/5u^ + 3/3t;^ F+. (4.35b) 

The compatibility conditions = F^ can also give rise to Eqn. (11.11) . Note that this form 
of Lax pair (14.341) has been presented in Ref . [l9j , and other three Lax pairs are given here 
for the first time. 
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To this stage, with the two-singular-manifold method, we have obtained the Lax pairs of 
the (2+l)-dimensional Gardner equation, i.e.. Systems (14.321) — ( [4. 35p . In the next section, 
we will construct the relationship between the singular manifolds 0, and eigenfunctions 



5. Relationship between the singular manifolds and eigenfunctions 

Using Eqns. fl4:8|) - fICTD and we rewrite Eqns. fICTD - flOTj) as 

(5.1) 
(5.2) 
(5.3) 
(5.4) 





4^xx 




i^x 


4>x 




i^tx 


_ fxx 


i^x 




^x 










rx 


4>x 






^x 


i^x 








= ? + 




^x 


(Px 






- ^ _ 


i^x (, 




y^x 


r V 



V'"^ ^ 0x ^x 0x ^x 

_2^ + ^ - ^ -3^ -3^^ (5 6) 

^~ <^x ^x 0x ^x) ' 

Integrating Eqns. (15.11) and (15. 2p with respect to x yields 

(/)^ = ^- '0+, (5.7) 

^x = ^tr- (5.8) 

Then, by substituting Eqns. (15. 7p and (15. 8p into Eqns. (15.30 — (15.60 . we obtain 

0j, = V^;V^+ + CV'x, (5.9) 

^y = i,^^r -rx^l. (5.10) 

<\>t = ^+ - 2 ^+ + 2 ^+ +4^; (5.11) 

= ^+ + 2 ^+ - 2 - 4 (5. 12) 

From Eqns. (15. 7p — ( 15.120 . it is shown that the singular manifolds and (p are determined 
by the eigenfunctions and ijj'^. So, by defining the singular manifolds and ip in the 
abbreviated form as jol 

= A(V'-,^+), (5.13) 

<^ = fi(^-,V'+), (5.14) 

Eqns. (15. 7p — (15.120 can be written as 
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[A(^-,^+)]^ = V'.-V^+, (5.15) 

[l](V^-,V^+)]^ = ^-^+, (5.16) 

[A(^-,V'+)]^ = + (5.17) 

[n {r, ^+)], = r - c i^t. (5.18) 

[A(^-, V^+)] ^ = ^+ ^- - 2 ^+ + 2 i^t ry + 4 C (5.19) 

[1] (^-, V^+)] ^ = r i^t + 2 ^+ - 2 ^+ - 4 C i^y- (5.20) 



Here, we can see that Eqns. (15. 15p — ( [5.201) establish the relationship between two singular 
manifolds 0, ip and eigenf unctions ?/'"'", and possess the following relation 

A(^/'-,V'+) + ^^(^",^+) =^"^+. (5.21) 



6. Binary Darboux transformation 

The Darboux transformation method is a powerful tool to get the analytical solutions 
for the integrable NPDEs 35|, l36|]. The most obvious advantage of this method lies in its 
iterative algorithm, which is purely algebraic and can be easily achieved on the symbolic 
computation system. By virtue of the SMM and Darboux transformation, starting from 
the seed solution and solving the corresponding linear equation or system, one can obtain 
wide classes of exact analytical solutions for a NPDE, such as the soliton solutions, periodic 



solutions and rational solutions 35, 



36 



37 



38 



4Q|. 



Based on the Lax pairs fl4.32p and (14.331) obtained in Section 3, we can construct the 
binary Darboux transformation of Eqn. (II .ip . Although Eqns. (12.40 and (12. 5p are considered 
as an auto-Backlund transformation, it is actually not convenient to generate more and more 
complicated solutions in a recursive manner, because it only involves the transformation for 
potentials. In comparison, the Darboux transformation not only has the potential transfor- 
mation, but also establishes the relationship between the new and old eigenfunctions. Next, 
we turn our attention to the transformation of eigenfunctions. 

Let us assume that 



/'- /- I •^l I •^2 

(j) if 

^ = ^1 + — + — , 

V 



(6.1) 
(6.2) 



where ipi and tpf respectively correspond to the solutions of Lax pairs (I4.32p and (I4.33p . while 
fi and Qi {i = 1,2) are two differentiable functions to be determined, the new eigenfunctions 
ip'~ and ijj'^ also satisfy Lax pairs (14.321) and (I4.33P except that {u,v) is replaced by (u',v'), 
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namely, 



^'v = -i'xx V'x , (6.3a) 



a 



( — ^ 6 Pu' -\ — a"^ u'"^ — 3 av' + 3 au'^ ijj^ 

\ 2 J 

+ 4V'l;x + ( 6tt^^' - ^ ) ^1;, (6.3b) 



and 



i^'y^ = '^'xl- [(^W -—] tp'^, (6.4a) 



a 



— b pu -\ — a u — 3av — 3au^ \ w~: 

2 J 

+ 4^1t.- f6a«'-^V:t (6.4b) 



Substituting Eqns. ^M)-UM and Ansatzs ([6lD-(l621) into Eqns. ([631) - and equating 
to zero the coefficients of like powers and (p, with symbolic computation, yield the following 
set of equations: 

hx = Bh - (6.5) 

/i. = |^/i + f (6-6) 

fu = ^h + ^(-^-AA^ + AA^-2aAu-AAr.-AAB)AQ.7) 
^ B \ Ax Ax a / 

h = 0, (6.8) 

and 

g2x = Ag2-'^^, (6.9) 

^72, = ^ + ^ {-i^ty + B i^tx) , (6.10) 

g.. = $g. + ^(-^+4B$-4B^- + 2aBu + 4Br,-4AB\i6.11) 
^1=0. (6.12) 

From Eqns. (16.51) — (16. 7p and (16. 9p — (16.111) . /i and g2 can be determined as 

/i = -V^-A(^r,V^+), (6.13) 
= (6.14) 

10 



Now, the relation between the new and old eigenfunctions has been constructed. Therefore, 
we arrive at the Darboux transformation for System (12. ip in the form 

2 



u 



u + 



v + 



a 
2 

a 



A(^-,^+)' 



jlog 


rA(V'-,V'+)" 


). 






jlog 


-A(^-,^+)j 






), 



(6.15) 
(6.16) 
(6.17) 
(6.18) 



7. Iteration of binary Darboux transformation and Grammian solutions 

In this section, we will perform the binary Darboux transformation N times and progres- 
sively generate the analytical Grammian solutions. After the second iteration of the binary 
Darboux transformation, the new potential functions and eigenfunctions are expressed as 



u[2] 
v[2\ 



with 



+ i 

a 

, 2 
v' ^- 
a 



m 
m 



u-\ — < log 

a 



V -\ — < log 

a 



m 
m 
m 
m 



A(r,^^ 



A(^'-,^'+), 



(7.1) 
(7.2) 



(7.3) 
(7.4) 



where 0' and ip' are the new singular manifolds for {u',v'), while ((t>,}p) and ($[2],E[2]) are 
the first and second iterative r-functions in the Hirota method 
and (17.41) . one can immediately calculate that 
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12|. FromEqns. ([731) 



m 



fi(v^r,^+ 



(7.5) 
(7.6) 



Taking and ipf' ( i = 1, 2, 3 ) as the solutions of the Lax pairs (14.321) and (I4.33P for {u, v) 
and iterating the Darboux transformation, we can get the following results: 



u[3] 
v[3] 



u + 



V + 



a 
2 
a 



-11 



(7.7) 
(7.8) 
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with 



m 



A(^3",v^2+; 



(7.9) 



(7.10) 



Following the same procedure above, we iterate the Darboux transformation times with 
symbolic computation and obtain 



u[N] 
v[N] 



u + 



V + 



a 
2 
a 



log 



(7.11) 
(7.12) 



with 



m] 



where ipi ip^ {i = 1, 2, ■ ■ 
In illustration, we take u 
yielding 



A(A-,^i+ 
A(V'2",V'i+ 

A(^^,^i+ 

fi(^r,^+ 



A(^r,^2+; 
A(^2",^2+; 



^(^2)^2 



A(^r,V'^) 
fi(^r,^+) 

^{'^2^i^N) 



(7.13) 



(7.14) 



, A^) satisfy Lax pairs (02]) and K3^ . 

= as the seed solutions for Lax pairs 04.321) and fl4.33p . 



^2" 



exp {ki X + hy + Wit} , 
exp {pi X + rriiy + riit} , 
exp {k2X + + W2t} , 
exp {p2 a; + m2 y + n2 1} , 



2/3 



f Pi + Pt 



Hi 



where = f A;, - A;2, = Akf - '-fkf + '-jfh, m, 
with fcj and pi {i = 1, 2 ) as arbitrary constants. 

Integration of Eqns. ( IS.lSp — ( l5.20p with respect to x, y and t results in 

h 



(7.15) 
(7.16) 
(7.17) 
(7.18) 

+ '-fpl + '-^P. 



h+Pi 



^1 i^i + ^1, 



(7.19) 
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k2 


+ P2 






"'I 


-\- 7)0 








+ Pi 






h 


+ Pl 




P2 




+ P2 




Pi 


h 


+ P2 




P2 


k2 


+ Pl 





+ S2, 


T 2 




.7/1+ 


+ O4, 


i't 










-53, 


i't 


- ^4, 



1,2,3,4) are all arbitrary constants. Particularly, setting 5i 



where 6i ( i 

^3 = 54 = 0, the r-functions can be expressed as follows: 



where 



F2 
R 



= l + A;iFi, 
ip = PiFi-l, 

$[2] = 1 + A;i Fi + k2 F2 + k^k2 R Fi F2, 
S[2] = 1 - Fi - p2 F2 + P1P2 R Fi F2, 



exp {{ki +pi)x + {h + nil) V + + ni) t + 61} , 

exp {{k2 +P2)x + {I2 + ^2) y + {w2 + n2) t + 62} , 
{ki - k2){pi -P2) 



(ki 

(fci+Pi)6» 



Pl)ik2+P2] 

= ki+pi, 



g-(fc2+P2)e2 _ ^ 



(7.20) 
(7.21) 

(7.22) 

(7.23) 
(7.24) 
(7.25) 
(7.26) 
(7.27) 

62 = I and 

(7.28) 
(7.29) 
(7.30) 
(7.31) 

(7.32) 
(7.33) 

(7.34) 

(7.35) 
(7.36) 



It is of interest to note that Expressions (17.281) — f l7.31l) are very similar to the forms of soliton 



solutions obtained with the Hirota method 



4l|. Substituting Expressions (17:281) - fl7:3T]) into 



sense of Refs. 





L6|) 


a 


41 



42 



]) and (17. 2p . one can get the one and two soliton solutions in the 



43| 



8. Conclusions 

The singular manifold method from the Painleve analysis plays a vital role in investigat- 
ing many important integrable properties for the NPDEs. In this paper, we have successfully 
applied the two-singular-manifold method to the (2+l)-dimensional Gardner equation and 
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derived its Hirota bilinear form, bilinear Backlund transformation, Lax pairs, as well as bi- 
nary Darboux transformation. With the help of symbolic computation, we have performed 
the A^-time iterative algorithm of binary Darboux transformation to generate the N x N 
Grammian solution. 
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Appendix: Identity properties of the bilinear operator 

The following identities are used in the derivation of the bilinear Backlund transforma- 
tion. 

Dx Co a(x) ■ a{x) = a{x) ■ cq a(x) = 0, ( cq = const. ), (A.l) 
{D^a-b)cd- ab{D^c- d) = {D^a ■ c)bd - ac{D^b ■ d) ^ D^{ad) ■ {be), (A.2) 
{Dla-b)cd-ab(Dlc-d) = (Dla ■ c)bd - ac{Dlb ■ d) - 2 D^ad ■ (D^b ■ c), (A.3) 
{Dla-b)cd- ab{Dlc- d) = {Dla ■ c)bd - ac{Dlb ■ d) -3D^{D^a- d) ■ (D^b-c), (A.4) 
{D^Dy a ■ b) cd — ab{ D^Dy c- d) — [D^Dy a-c)bd — ac{ D^Dy b • d) 

-D^{ad) ■ (Dyb-c) - Dy{ad) -{D^b-c). (A.5) 

References 

[1] J. Weiss, M. Tabor and G. Garnevale, J. Math. Phys. 24, 522, 1983. 

[2] J. Weiss, J. Math. Phys. 24, 1405, 1983. 

[3] M. Musette and R. Gonte, J. Math. Phys. 32, 1450, 1991. 

[4] A. G. NeweU, M. Tabor and Y. B. Zeng, Physica D 29, 1, 1987. 

[5] P. G. Estevez, Phys. Lett. A 171, 259, 1992. 

[6] W. H. Steeb and N. Euler, Nonlinear Evolution Equation and Painleve Test (World 
Scientific, Hong Kong, 1988). 

[7] F. Gariello and M. Tabor, Physica D 39, 77, 1989. 

[8] M. Musette and R. Gonte, J. Phys. A 27, 3895, 1994. 



14 



[9] P. G. Estevez and P. R. Gordoa, Inverse Problems 13, 939, 1997. 
[10] P. G. Estevez, P. R. Gordoa, L. M. Alonso and E. M. Reus, J. Phys. A 26, 1915, 1993. 
[11] R. Conte, M. Musette and A. Pickering, J. Phys. A 28, 179, 1995. 
[12] P. G. Estevez and P. R. Gordoa, Teor. Matem. Fizika 99, 370, 1994. 
[13] J. M. Cervero and P. G. Estevez, J. Math. Phys. 39, 2800, 1998. 
[14] P. G. Estevez, J. Math. Phys. 40, 1406, 1999. 
[15] G. Das and J. Sarma, Phys. Plasmas 6, 4394, 1999. 

[16] M. P. Barnett, J. F. Capitani, J. Von Zur Gathen and J. Gerhard, Int. J. Quantum 
Chem. 100, 80, 2004. 

[17] B. Tian and Y. T. Gao, Phys. Lett. A 342, 228, 2005. 

[18] W. P. Hong, Phys. Lett. A 361, 520, 2007; Y. T. Gao and B. Tian, Phys. Plasmas 13, 
112901, 2006; Phys. Plasmas (Lett.) 13, 120703, 2006; Phys. Lett. A 349, 314, 2006; Y. 
T. Gao, B. Tian and C. Y. Zhang, Acta Mech. 182, 17, 2006. 

[19] B. G. Konopelchenko and V. G. Dubrovsky, Phys. Lett. A 102, 15, 1984. 

[20] B. G. Konopelchenko, Inverse Problems 7, 739, 1991. 

[21] X. G. Geng and C. W. Gao, Nonlinearity 14, 1433, 2001. 

[22] I. Anders, Asymptotic Analysis 19, 185, 1999. 

[23] Y. Chen and Z. Y. Yan, Chaos, Solitons and Fractals 26, 399, 2005. 
[24] G. F. Yu and H. W. Tarn, J. Math. Anal. Appl. 330, 989, 2007. 

[25] R. Hirota, The Direct Method in Soliton Theory (Cambridge Univ. Press, Cambridge, 
2004). 

[26] A. Nakamura, J. Math. Phys. 22, 1608, 1981. 

[27] C. Rogers and W. Shadwick, Bdcklund Transformations and Their Applications (Aca- 
demic Press, New York, 1982). 

[28] V. E. Zakharov, What is Integrability (Springer Press, Berhn, 1991). 

[29] C. Y. Zhang, Y. T. Gao, X. H. Meng, J. Li, T. Xu, G. M. Wei and H. W. Zhu, J. Phys. 
A 39, 14353, 2006. 

15 



[30] H. Q. Zhang, X. H. Meng, T. Xu, L. L. Li and B. Tian, Phys. Scr. 75, 537, 2007. 

[31] B. Tian, Y. T. Gao and H. W. Zhu, Phys. Lett. A 366, 223, 2007. 

[32] S. Ganesan and M. Lakshmanan, J. Phys. A 20, L1143, 1987. 

[33] R. Radhakrishnan, M. Lakshmanan and M. Daniel, J. Phys. A 28, 7299, 1995. 

[34] T. Alagesan, Y. Chung and K. Nakkeeran, Chaos, Solitons and Fractals 26, 1203, 2005. 

[35] V. B. Matveev and M. A. Salle, Darboux Transformations and Solitons (Springer Press, 
Berlin, 1991). 

[36] C. H. Gu, H. S. Hu and Z. X. Zhou, Darboux Transformation in Soliton Theory and its 
Geometric Applications (Shanghai Scientific and Technical Publishers, Shanghai, 2005). 

[37] Y. S. Li and J. E. Zhang, Chaos, Solitons and Fractals 16, 271, 2003. 

[38] Z. X. Zhou, J. Math. Phys. 39, 986, 1998. 

[39] Q. H. Park and H. J. Shin, Physica D 157, 1, 2001. 

[40] M. J. Ablowitz, S. Chakravarty, A. D. Trubatch and J. Villarroel, Phys. Lett. A 267, 
132, 2000. 

[41] J. Hietarinta, J. Math. Phys. 28, 2094, 1987. 
[42] S. F. Deng, J. Phys. A 39, 14929, 2006. 
[43] S. F. Deng, Phys. Lett. A 362, 198, 2007. 



16 



